Criterion for bosonic superfluidity in an optical lattice 
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We show that the current method of determining superfluidity in optical lattices based on a visibly 
sharp bosonic momentum distribution n(k) can be misleading, for even a normal Bose gas can have 
a similarly sharp n(k). We show that superfluidity in a homogeneous system can be detected from 
the so-called visibility (v) of n(k) — that v must be 1 within 0{N~^^^), where A*' is the number of 
bosons. We also show that the T = visibility of trapped lattice bosons is far higher than what is 
obtained in some current experiments, suggesting strong temperature effects and that these states 
can be normal. These normal states allow one to explore the physics in the quantum critical region. 



There has been strong interest in using cold atoms 
in optical lattices to simulate strongly correlated many- 
body systems so as to shed light on many long stand- 
ing problems in condensed matter physics. The in- 
terest began a few years ago with experiments on the 
superfluid-insulator transition of bosonic atoms in opti- 
cal lattices fl|, and has grown rapidly since the achieve- 
ment of fermion pair condensation near a Feshbach res- 
onance 0. One class of very important problems, in- 
cluding high Tc superconductivity, is understanding how 
superfluid order (bosonic or fermionic) develops, and how 
the superfluid transforms into other correlated many- 
body states as the interaction parameters are varied. To 
achieve this goal, it is necessary to reach quantum degen- 
eracy in a lattice, and to identify the presence of super- 
fluidity. 

At present, the method commonly used for identifying 
superfluidity of bosons is through the "sharpness" of the 
diffraction spots in their momentum distribution n(k) 
P, 0, H, [^. Despite its popularity, there has been no 
effort to characterize this "sharpness" precisely. As far 
as we can tell, a peak is considered "sharp" if its width is 
visually much smaller than the separation of peaks. To 
be sure, a macroscopic bosonic superfluid is characterized 
by a (5-function peak in n(k) (of order N where N is the 
number of particles). Unfortunately, the presence of such 
(5-function is hard to discover due to finite experimental 
resolution. Instead, one relies on the estimate of sharp- 
ness mentioned above, which is consistent with but not 
a proof of superfluid correlation, as we explain below. 

The purpose of this paper is to point out a number of 
facts crucial for identifying superfluid order for bosons in 
optical lattices. We show that (I) even a normal Bose gas 
above Tc can have a diffraction pattern as sharp as those 
in current experiments. Identifying superfluid order from 
the sharpness of n(k) as practiced today is therefore un- 
reliable. (II) For homogeneous systems, the presence of 
superfluid order implies that the so-called "visibility" (w) 
must be 1 within 0{N^^/^). We also present (III) the 
visibility at T = as a function of lattice parameters for 
the "wedding cake" structure of harmonically confined 
lattice bosons. In this case, v deviates from 1 when the 
superfluid regions are sufficiently small. In current exper- 



iments, this typically occurs after more than one Mott 
layers have developed. Because of the high sensitivity 
of V to superfluid order, this visibility curve is a good 
calibration of temperature effects in the system. These 
results have strong implications for the interpretation of 
many current experiments, discussed at the end. 

In our discussions, we shall use the identification 
adopted in all current experiments P, [3, H, [B] that the 
observed diffraction pattern is related to the momentum 
distribution of the system through a ballistic expansion 
of the cloud (see eq. ^ below) . While this has not been 
proven rigorously, it is consistent with the fact that the 
confinement energy of a Wannier state in the tight bind- 
ing limit is much larger than the interaction energy [2] . 

(A) Normal Bose gas in a lattice: Consider an 
ideal Bose gas with N bosons in an optical lattice with 
volume ri. The Hamiltonian is H = * = 

h, = ^{h^/2m)df + V^sin^{TTXi/d), Va > 0. When 
Vo is sufficiently large, only the lowest band (with en- 
ergies i?k = ~2t ^-^^ ^ ^ cos fcid and Bloch functions 
^'k(x)) is thermally occupied. Above the Bose conden- 
sation temperature the chemical potential is de- 
termined by N/fl = X^k /_B('^k), where fg{x) = 
(e(2:-M)/fcBT _ i)-i^ and Y^'^ is a sum over the first Bril- 
louin Zone. At Tc, fi reaches the bottom of the band 
Eq. The momentum distribution for T > Tc is ri(q) = 
X^k /s(^k)|*I'k(q)P, where ^k(q) is the Fourier trans- 
form of ^'k(x). Experimentally, one measures the column 
distribution, Afj^(qj_) — J dg2^^(q), where q = {qj^,q^). 
Since 5'k(q) is non-zero only when q = k -I- G, where G 
is a reciprocal lattice vector, and since Ey^ = T'k+c 
have N^{qj_) = /dg,/B(Sq)|*q_G(q)|', with q - G in 
the first Brillouin zone. For a narrow band, |^q_G(Q)P 
is a Gaussian centered at q = decaying on the scale 
27r/d. Hence iVj^(qj^) is composed of peaks centered at 
reciprocal lattice vectors G, with the shape of the peak 
given entirely by the variation of /g(_Bq_G) around G 
and an overall envelope given by |^'q_G('3)P- 

In Fig. 1(A) we show iVj_(q_L) for a lattice with one 
boson per site {N/fl = d~^) and Vo — l5Eji at a tem- 
perature T = l.lTc, where En = %^ tt"^ / {2m(P) is a "re- 
coil" energy. We have found numerically that in this 
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FIG. 1: (A): A^^(q) for an ideal Bose gas with one particle 
per site in an optical lattice at T = 1.1 Tc for Vo = IbEn. A 
and B refer to those vectors defined in Eq. ([ij. (B): A'^^(q) 
along the g^^-axis. (C): v vs T/Tc for this ideal Bose gas. 



system ksTc = 0.455, where B ~ 12t is the bandwidth 
of the lowest band, thus fcsT < B. As we shall see, 
this leads to sharp peaks distributed on a 2D square lat- 
tice with spacing 2TT/d. Fig. 1(B) shows that the peaks 
are visibly "sharp" , with a full width at half maximum 
{Aq)T = 0.1(27r/(i). This demonstrates that diffraction 
spots in n(q) with width much less than a reciprocal lat- 
tice spacing is not proof of Bose condensation. 

(B) Condition for quantum degeneracy: The re- 
sults in (A) can be understood by considering the con- 
dition for quantum degeneracy. When k^T < B, the 
most thermally occupied states are near the bottom of 
the energy band, for which we can use the approximate 



spectrum Ey^ sa ~6t + h k /2m* , where m* is the effec- 
tive mass defined as (to* = h^/{2td'^)). The "lattice" 
thermal wavelength Ay — h / ^/2Trm*kBT is reduced from 
the free space value A^' by ^rajm* . The condition for 
Bose condensation for one boson per site, which is also 
the condition for quantum degeneracy (At ~ , becomes 
ksTc = 0.55i?. The difference from the numerical result 
(0.45i?) is due to the effective mass approximation. The 
width of the spot is in general proportional to A^^. 

Note that the change of thermal wavelength means 
that in a lattice is reduced by a factor of m/m* . This 
poses a severe challenge to reaching quantum degener- 
acy in the deep lattice limit. For V^/Ej^ — 10, 15, 30, we 
have m/m* = 0.25, 0.09, 0.007 respectively. Without lat- 
tice, for gases with 10^ bosons in harmonic traps, Tc is 
typically 10~^K and the lowest temperature attainable 
today is 10~^K. For deep lattices with m/m* ^ 10~^, 
one can barely reach quantum degeneracy even at the 
lowest temperature attainable today 8]. 

(C) Visibility and Bose condensation: The 
visibility, originally introduced to study short range 
coherence 0, is defined as 



Na-Nb 



(1) 



where Na = NJGx), Nb = N^{Gh), G = 27r/d, Gx is a 
reciprocal lattice vector; Gh is Gx rotated by 45° around 
the z axis and is not a reciprocal lattice vector. In the 
superfluid phase, Nj_{Gx) N while Nj^{Gn) ^ N^^^ 
(see later discussions), so we have v ^ 1. The visibility 
of an ideal Bose gas in a lattice with Vo — 15i?fl and 
one boson per site is shown in Fig. 1(C). The visibility 
is 100% at r < Tc but decreases sharply above Tc. It is 
interesting to note that despite its sharp drop at Tc, v 
decays slowly, remaining at 0.1 at T = lOTc. 

For interacting bosons in a sufhciently deep lattice, 
the bosons are confined to the lowest band, described 
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R"R' ~r 

means near- 



by the Bose Hubbard model H - 

h.c.) + ^X]r'*r("-r ~ 1)' where (R, R', 
est neighbors, creates a boson in the Wannier state 
WR(r) = T-3/2^^g-»k R^j^(-j.) located at site R, T^ is 

the number of lattice sites, and riR = a^a^^. The hop- 
ping integral t and the interaction parameter U are calcu- 
lated from the eigenstates of hi and the s-wave scattering 
length in a straightforward manner. Since the Fourier 



transform of WR(r) is of the form Wj^(q) 
the momentum distribution is 



-iq-R 



.(q) = |u;(q)P ^(at,ai,,)e^'i-(^-^') 



w(q). 



(2) 



R.R' 



In a homogeneous superfluid, {aT^a-^,) is essentially given 
by l^fp for R 9^ R', * = (orHSJ- Denoting the number 
of condensed bosons as Nn = i'^l^'P Il0|, we have 



.(q)=[(A^-7Vo) + |vI/p/(q)] Wq)!^ 



(3) 



where /(q) — \ X]rS"*^P- For a narrow band, w(q) is 
well approximated by |w(q)P = 11^=2;.^,^ W((ji), >V(fc) 
g-fc /a fj ^ ijfi^ cubic lattice, we have 

/(q) = Y{^F{q^). F{k) = [sin(TA:d/2)/sin(fcd/2)]2, which 
peaks sharply at reciprocal lattice vectors G with a width 
'-^ n/Ld. Since the product yV((7a;)W('Zy) has the same 
value at Gx and Gn and since T(0)T(27r/c?) — L*, it is 
simple to show from Eqs. (H]) and ^ that 



1*1 



2 / W{q,)dq, 



|*|2 + (iV-iV„)y' JWiq,)F{q,)dq,L* 



(4) 



Simple integration shows that y ~ ^T^- Since j'l'p = 
No/L^, we then have in the superfluid phase., u = 1 — 
0(1/A^2/^) nl. In order for the visibility to deviate from 
1 by a non-zero but small amount, f = 1 — e, we need 
No/N - l/(eT2) ^ l/(e7V2/3)^ Hence, even if the con- 
densate fraction is very small, as long as it is larger than 
l/(eAf2/^), the visibility is essentially l,(l>w>l — eto 
be precise). For example, a visibility of 0.95 for a system 
with 10^ bosons with about a few bosons per site means 
No/N « 10-3. 

(D) Visibility of lattice bosons in a harmonic 
trap: Many recent experiments investigating quantum 




FIG. 2: (A): Schematic phase diagram for a homogenous lattice Bose gas [T^. The superfluid, (quantum critical) normal, and 
Mott insulator phases are labeled SF, N, and M respectively. The solid blue line is the critical temperature for the system. 
The shaded area is a crossover. I and II represent different physical processes. (B) : Density vs radial distance for lattice bosons 
in a harmonic trap for conditions in ref. at Vq/Er = 12 (dotted black curve), 14 (dot-dashed red), 16 (dashed green), and 
20 (solid brown). Inset: phase diagram for a homogeneous system. The region outside the lobes is superfluid. Vertical lines 
represent values taken by fj,{r) for the corresponding Vo- (C): T — visibility vs Vq. For Vb < 14.7i5fl, where the second Mott 
shell begins to appear, the system has a large superfluid core, which yields v = 1. The upper solid red (lower, solid blue) curve 
is the visibility when the superfluids separated by the Mott shell are in (out of) phase. The circles are experimental data from 
ref . [l3| . The dashed line is a calculation including short range coherence assuming all superfluid regions have become Mott 
phases. The dotted line is the superfluid fraction at T = 0. 



phase transitions of a lattice Bose gas or Bose- Fermi mix- 
tures have measured the visibiHty of these systems in har- 
monic traps as a function of lattice height 14 0, 3, d, [l^ ■ 
Except for a single case in ref. d {Vo = 5Er) which finds 
V = 1 in the regime where a majority of bosons should 
be superfluid in the ground state, the data reported in 
ref . [1, 0, H, shows that v < 0.8 in a similar regime. If 
there was no harmonic trap, a visibility ?; = 0.8 means 
the lattice Bose gas is normal, as shown in Section (C). 
On the other hand, in a harmonic potential V(r), it is 
well known that the system develops alternating layers of 
superfluid and Mott phases (the so-called "wedding cake" 
structure). When a sufficiently large number of bosons is 
converted from the superfluid phase to the Mott phase, 
the visibility will begin to drop. In addition, finite tem- 
perature or heating effects can also destroy phase coher- 
ence and reduce visibility. 

To understand the general behavior of the visibility, let 
us consider a region in the harmonic trap (say, around 
r) where the lattice Bose gas turns into a Mott phase as 
the lattice depth Vo (and hence the ratio U/t) increases. 
Within the local density approximation (LDA), we can 
treat this region as a bulk system for which the physical 
process is represented as a path in the phase diagram of a 
homogenous lattice gas shown schematically in Fig. 2(A), 
which plots the transition temperature Tc as a function 
of U/t. At U/t = 0, Tc is given by the quantum degener- 
acy condition discussed in Section (B). It drops to zero 
at the quantum critical point {U/t)*. The shaded line 
in Fig. 2(A) is the crossover from the quantum critical 
region (a normal phase with no clear sign of a gap) to the 
Mott region (a normal phase with an interaction gap). 

Since experiments are performed at finite temperature, 
any physical trajectory connecting the superfluid phase 



to the Mott phase must pass through the quantum crit- 
ical region. Typically, as K increases, the system heats 
up due to a variety of reasons : spontaneous emission, 
tiny vibrations of the apparatus, etc. The physical pro- 
cesses may therefore look like trajectories / or // shown 
in Fig. 2(A). The states (a) and (6) are in the super- 
fluid phase. The final states (c) and (e) are in the Mott 
regime. The state {d) is in the normal regime. The prox- 
imity to a quantum phase transition can be measured by 
the length of the trajectory passing through the quan- 
tum critical region. For example, process / is close to 
the quantum phase transition, whereas // is not. In ho- 
mogeneous systems, for both I and II one starts off with 

V = 1 and a sharp momentum distribution n(k) in the su- 
perfluid region. For /, v drops sharply and n(k) becomes 
blurry quickly across the transition point {U/t)* . For //, 

V drops slowly as the system leaves the superfluid phase, 
and n(k) remains sharp in the quantum critical regime 
close to Tj,. 

The proximity to a quantum phase transition can also 
be estimated by comparing the measured visibility to 
the T = visibility calculated using standard mean 
field methodlil and LDA. To be concrete, we focus on 
the system in ref.[l2| because it has the most detailed 
analysis of data among current experiments on lattice 
bosons 14i]. The physics illustrated here, however, should 
be applicable to boson- fermion mixtures 0, S]. We begin 
by calculating the order parameter {a-^ = ^'r and den- 
sity {n-^ of an infinite lattice as a function of chemical 
potential /i and the interaction ratio t/U for a homo- 
geneous system. The phase boundary between super- 
fluid and Mott phases is a sequence of "Mott lobes" as 
shown in the inset of Fig. 2(B). The regions within differ- 
ent lobes are Mott phases with different (integer) num- 
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ber of bosons per site. In a trap V{r), both and 
(n-p^) are position dependent, since becomes (within 
LDA) /i(r) = /i — y(r). In this way, we obtain the 
density profiles in Fig. 2(B). (We mention that our den- 
sity profiles differ from those in ref. [l^l, see [3]. These 
differences, however, will not affect our points below.) 
We then calculate v from Eqs. ([1]), ([2]) with (a^Oj^,) = 
('^r) R' + ^R^R' ( 1 ^ '^R R' ) ■ If R' are in discon- 

nected superfluid regions, the product ^'r^'r/ depends 
on the relative phase A9 between these regions. 

In Fig. 2(C), we have plotted our result for i; as a 
function of Vo for the system in ref.[l3]- The differ- 
ent curves correspond to different ways of treating the 
relative phase A9 between different disconnected super- 
fluid regions. One sees that the T = visibility differs 
strongly from the experimental data (shown as circles in 
the figure). This large difference, however, is not due to 
the differences in our wedding cake structures. This is 
because both structures contain a superfluid core below 
Vo = 14.7i?jj (where the second Mott shell begins to de- 
velop) that is so large that v must be 1 as long as the su- 
perfluid is not destroyed. The disagreement with experi- 
ments implies that all superfluid regions that should exist 
atT = have turned normal (hence the much weaker vis- 
ibility) , which can only occur if the temperature is above 
Tc in these regions (i.e. the system that should be in the 
superfluid state (6) in Fig. 2(A) at T = is found to be 
in state (d) above Tc). The physical process is therefore 
quite far from the quantum critical trajectory [3]. 

In order to account for the visibility deep in the Mott 
regime, the authors of ref.[l3| considered short range cor- 
relations in a perturbative manner and found good agree- 
ment with their data, provided one makes the assump- 
tion that all the superfluid regions are converted into 
the Mott phase. We have repeated this procedure with 
our wedding cake structure [l5| and have obtained sim- 
ilar agreement (dashed line in Fig. 2(C)). The assump- 
tion that lead to this agreement, which eliminates all 
contributions from superfluid to visibility, is consistent 
with the picture that all superfluid regions has gone nor- 
mal due to temperature effects. We would like to point 
out that our mean field calculations do not include these 
short range coherence but that their inclusion would only 
raise the visibility curves in Fig. 2(C) to even larger 
values. Finally, it is instructive to look at the conden- 
sate fraction Nq/N at T = a,s a. function of Vb, where 
No = Er I*rP- We see from Fig. 2(C) that at T = a 
visibility as high as 0.8 (which occurs when Vo > U.TEn) 
represents a condensate fraction No/N ^ 0.05. 

(E) Implications for recent experiments: Just as 
in ref.[l3, ref.[l,[l] also show visibility v ~ 0.8 for lattice 
heights where the system should be superfluid at T = 0. 
The physical processes in ref. 0, Q are therefore quite far 
from being a quantum phase transition (QPT). (To show 
that a physical trajectory is close to a QPT, it is neces- 
sary to demonstrate that the quantum critical region tra- 



versed in the process is very narrow). Moreover, the de- 
crease of visibility of the bosons when fermions are added 
suggests that fermions may be increasing the tempera- 
ture of bosons. Our discussions in Section (A) and (B) 
also show that the recent claim of observation of super- 
fluid correlation of fermions in an optical lattice Q based 
on the sharpness of n(k) is not conclusive. The claim 
would have been established if the bosonic molecules af- 
ter the sweep were found to have visibility v — 1. 

Our study indicates that the problem of heating is 
prevalent in current experiments. We hope our findings 
will stimulate serious efforts to determine the tempera- 
ture of lattice gases and more rigorous ways to achieve 
quantum degeneracy in lattices. This work is supported 
by NSF Grant DMR-0426149 and PHY-0555576. 
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